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We report measurements of the third harmonic coefficient of the azimuthal anisotropy, V3, known 
as triangular flow. The analysis is for charged particles in Au+Au collisions at ^snn = 200 GeV, 
based on data from the STAR experiment at the Relativistic Heavy Ion Collider. Two-particle 
correlations as a function of their pseudorapidity separation are fit with narrow and wide Gaussians. 
Measurements of triangular flow are extracted from the wide Gaussian, from two-particle cumulants 
with a pseudorapidity gap, and also from event plane analysis methods with a large pseudorapidity 
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gap between the particles and the event plane. These results are reported as a function of transverse 
momentum and centrality. A large dependence on the pseudorapidity gap is found. Results are 
compared with other experiments and model calculations. 

PACS numbers: 25.75.Ld 



I. INTRODUCTION 

The study of azimuthal anisotropy, based on Fourier 
coefficients, is recognized as an important tool to probe 
the hot, dense matter created in heavy-ion collisions [U 
2 . The first harmonic coefficient v\ , called directed flow, 
and the second harmonic coefficient V2, called elliptic 
flow, have been extensively studied both experimentally 
and theoretically, while higher even-order harmonics have 
also garnered some attention [3]. In contrast, odd har- 
monics of order three and above were overlooked until 
recently. This is because in a picture with smooth ini- 
tial overlap geometry, it had been assumed that higher- 
order odd harmonics are required to be zero by symme- 
try. More recently it has been realized that event-by- 
event fluctuations break this symmetry [4]-[6]. The event 
plane approximates the participant plane and for rea- 
sonable event-plane resolutions the measured v n are not 
the mean values, but closer to the root-mean-square val- 
ues [7]. As a consequence, higher-order odd harmonics 
carry valuable information about "hot-spots" or "lumpi- 
ness" in the initial state of the colliding system [BlU6|. 
The third harmonic coefficient - sometimes called trian- 
gular flow, but probably not related to triangular con- 
figurations in the initial state - is thus a new tool to 
study initial state fluctuations, the subsequent evolution 
of the collision system, and possibly the production of the 
near-side ridge [H HZ] • Theoretical studies suggest that 
«3 is more sensitive to viscous effects than vi because 
finer details of the higher harmonics are smoothed more 
by viscosity [TUj. It also appears that the mean value 
of the triangular eccentricity from central to mid-central 
collisions is independent of the geometric model used for 
the initial overlap [TH] , unlike elliptic eccentricity. This is 
probably because V3 is an odd harmonic and dominated 
by fluctuations. Rapidity-even v\ is also dominated by 
fluctuations, but is complicated by the correction needed 
for conservation of momentum |19| . Higher odd harmon- 
ics are thought to be less useful because of non-linear 
terms coming from lower harmonics |20) . Thus V3 is an 
ideal flow harmonic to study viscosity because it is al- 
most insensitive to the model for the initial conditions 
and more sensitive to viscosity. 

In order to separate the long-range correlations 
of interest from short-range correlations, we present 
measurements, based on the azimuthal angle 0, of 
(cos[3(0j — 4>i)])i^j vs. the pseudorapidity separation 
A77 = rji — rjj between the two particles, fit with narrow 
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FIG. 1: (Color online) «|{2} vs. transverse momentum and 
pseudorapidity separation for charged hadron pairs in 200 
GeV Au+Au minimum-bias collisions. The pr is for one of 
the particles, integrated over the pr values of the other par- 
ticle. 

and wide Gaussians. We present results derived from the 
wide Gaussian, for two-particle cumulants, and for the 
standard event plane methods, as a function of transverse 
momentum pt, pseudorapidity gap A77 , and centrality. 
An important experimental variable is the pseudorapid- 
ity gap between the particles being correlated. We com- 
pare our results to other experiments, and both transport 
and hydrodynamic models. 



II. EXPERIMENT 

About ten million Au+Au collisions at ^/s^n = 200 
GeV have been used in this study, all acquired in the 
year 2004 using the STAR detector with a minimum-bias 
trigger. The main Time Projection Chamber (TPC) [2T] 
of STAR covers pseudorapidity 1 77 1 <; 1.0, while two For- 
ward Time Projection Chambers (FTPCs) [H] cover 
2.5<|?7|<4.0. The extended range in r\ of the FTPCs was 
important because the analyses were done as a function 
of the 77 gap between particles. This requirement limited 
the study to the data collection years when the FTPCs 
were operational. The centrality definition of an event 
is based on the number of charged tracks in the TPC 
with track quality cuts of |r?|<0.5, a distance of closest 
approach (DCA) to the primary vertex less than 3 cm, 
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and 15 or more space points out of a total of 45. This 
analysis used events with vertex z coordinate (along the 
beam direction) within 30 cm from the center of the TPC. 
For each centrality bin, the number of participants and 
binary collisions can be found in Table III of Ref . [53] . 



III. ANALYSIS METHODS 

A. Event Planes 

In the standard event plane method 24J for u 3 , we re- 
construct a third harmonic event plane from TPC tracks 
and also from FTPC tracks. For event plane reconstruc- 
tion, we use tracks with transverse momentum p<r>0.15 
GeV/c, that pass within 3 cm of the primary vertex, 
and have at least 15 space points in the TPC accep- 
tance ( 1 77 1 < 1 . 0) or 5 space points in the FTPC accep- 
tance (2.5<|7j|<4.0). It is also required that the ratio 
of the number of actual space points to the maximum 
possible number of space points along each track's tra- 
jectory be greater than 0.52. In event plane calculations, 
tracks have a weighting factor w = pt in units of GeV/c 
for p T <2 GeV/c, and w = 2 GeV/c for p T > 2 GcV/c. 
Although the STAR detector has good azimuthal sym- 
metry, small acceptance effects in the calculation of the 
event plane azimuth were removed by the method of shift- 
ing pp. When using the TPC event plane, we used the 77 
sub-event method which provides an 77 gap, with an addi- 
tional small T) gap of ±0.05 between the sub-events [24] . 
The 77 sub-event method avoids self-correlations because 
the particles and the event plane are in opposite hemi- 
spheres. When using the FTPCs, we obtained the sub- 
event plane resolution from the correlation of the two 
FTPCs, but then used the full event plane from both 
FTPCs [23]. This introduced a large r\ gap between the 
particles in the TPC and the FTPC event planes. Since 
there is no overlap between the coverage of the TPC and 
FTPCs, there is no possibility of self-correlation when 
using the FTPC event plane. 



B. 2-particle correlations 

We studied v\{2} = (cos[3(0j — <f>i)])i^j vs. 
A77 between the two particles. For this two-particle 
cumulant method acceptance correction terms, 

which were generally small, were evaluated and applied. 
Figure [T] shows that there is a sharp peak for tracks 
close in r\ and at low px- This has also been seen by 
PHOBOS [M]. Our distribution of uf{2} vs. A77 can 
be well described by wide and narrow Gaussian peaks 
as shown in Fig. [2] for two centrality intervals. Using 
two Gaussians plus a flat background gave the same 
results for v% when integrated for all accepted pairs 
within the range |A?7|<2, as described below. The 
narrow Gaussian is identified as short range non-flow 
correlations like the Bose-Einstein correlation, resonance 
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FIG. 2: (Color online) «|{2} vs. the pseudorapidity sep- 
aration of the particles in pairs for charged hadrons with 
0.15<pt<2.0 GeV/c within two centrality intervals in 200 
GeV Au+Au collisions. Data are fit with narrow and wide 
Gaussians. Like Sign (LS), Unlike Sign (US), and Charge In- 
dependent (CI) cases are shown with only statistical errors. 
The dashed curves under the peaks are the wide Gaussian 
fits. 



decay, and Coulomb interactions, reduced by effects 
from track merging. The narrow peak disappears above 
Pt>0-8 GeV/c, so is unlikely to be from jet correlations. 
The wide Gaussian is the signal of interest in this 
paper and its fit parameters are used to calculate w|{2} 
as a function of centrality and transverse momentum 
for accepted pairs within the range |Afy|<2. For v 3 
integrated over pt and 77, 



<-3 2 {2}> = 



J b a vmWd(A v ) 
J b a Wd(A v ) 



(1) 



where W equals dN/ d(Arj) when weighted with the num- 
ber of particle pairs. The integration ranges for numera- 
tor and denominator are the same. To evaluate the effect 
of weighting we also used unit weights. The differential 
v?,{2}{pt) can be obtained from the scalar product [T] 
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FIG. 3: (Color online) The width and amplitude of the charge 
independent wide Gaussian as a function of transverse mo- 
mentum for most central (0-5%) and mid-central (30-40%) 
Au+Au collisions at t/snn = 200 GeV. The plotted errors are 
statistical. 
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v 3 {2}( PT ) 



(cos[3(0,-(pr) - (/),)}) 



(2) 



where the j th particle is selected from the px bin of in- 
terest. 

Figure [3] shows the pt dependence of the width and 
amplitude of the wide Gaussian fit to the data in Fig. [2] 
without a constant offset. Shown are results for the 0-5% 
most central and 30-40% mid-central collisions. Above 
0.8 GeV/c the distribution can be described by a single 
wide Gaussian. The amplitude increases with px and 
then saturates around 3 GeV/c. The px dependence of 
the width depends on centrality, with the 0-5% most 
central data showing first an increase in the width and 
then a gradual decrease, while for the 30-40% central 
data the width appears to gradually decrease for all px- 

Figure [4] shows the centrality dependence of the width 
and amplitude of the wide Gaussian. In peripheral colli- 
sions, the Gaussian width is narrow and well constrained 
by the data. As the collisions become more central, the 
width broadens reaching beyond 1.5 units in pseudora- 
pidity in the centrality range 10-40%. When the width of 
the wide Gaussian becomes broader than Arj = 1, it be- 
comes difficult, with the data from the TPC alone, to dis- 
tinguish between various functional forms for V3{2}(Ar/). 
The data points in the figure show the results when fitting 
a single wide Gaussian to the TPC data alone. The upper 
edge of the systematic error band corresponds to a width 
that would allow the fit function to extend out far enough 
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FIG. 4: (Color online) The width and amplitude of the wide 
Gaussian as a function of centrality for Charge Independent 
and Like Sign particles with 0.15<Pt<2.0 GeV/c for Au+Au 
collisions at ^/snn = 200 GeV. The errors on the data points 
are statistical. The upper edge of the systematic error band 
shows the width of the wide Gaussian required to also fit the 
data from the FTPC. 



to match the FTPC data at (Arj) = 3.21. On the other 
hand, if we include a constant background, we can also 
match the FTPC and TPC data with a wide Gaussian 
width consistent with the lower edge of the error band 
in Fig. |4j A larger acceptance in r\ is required to better 
constrain the functional form. Such a constraint could 
help distinguish between different physical mechanisms 
underlying the signal, such as stochastic fluctuations in 
the hydrodynamic phase |27j or decoherence of flux-tube 
like structures in the longitudinal direction [28j . 

Whether using the TPC data only or also including 
the FTPC data, the width of the wide Gaussian peak 
tends to become narrow in the most central collisions. 
The rise and then fall of the width of t>3{2}(A?7) mimics 
the rise and fall of the low px ridge amplitude reported 
in Ref. Reference [5] describes this centrality trend 
in terms of participant eccentricity fluctuations, where 
the fluctuations in mid-central collisions are well above 
statistical expectations. This can be attributed to the 
asymmetry of the overlap region which allows a nucleon 
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FIG. 5: (Color online) The third harmonic coefficient as a 
function of pseudorapidity for different centralities for Au+Au 
collisions at ^/s/vjv = 200 GeV, with track selection in the 
TPC of 0.15<p T <2.0 GeV/c. Results are shown for the event 
plane constructed either in the TPC or in the FTPCs. The 
horizontal lines are fits to the FTPC results. 



on the periphery of one nucleus to impinge on many nu- 
cleons in the center of the other nucleus thus amplifying 
the effect of fluctuations of nucleon positions in the pe- 
riphery of the nucleus. Thus the width of U3{2}(A?7) and 
the amplitude of the low pr ridge may be related to the 
same fluctuations. 
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FIG. 6: (Color online) The top panels show third harmonic 
coefficient as a function of pr for the wide Gaussian method 
and for the event plane in the TPC, for two centralities for 
Au+Au collisions at ^/snn = 200 GeV, for tracks in the TPC 
with — 1.0<77<1.0. The wide Gaussian was weighted with 
either the number of particle pairs or by unity. The bottom 
panels show the ratio of V3 from the wide Gaussian method 
to i>3 from the TPC sub-event method. 



With the event plane in the FTPCs there is a large 77 
gap between the particles and the plane, and V3 is flat 
for all centralities. This flatness means that acceptance 
effects at the edges of the TPC are not significant. Thus, 
even though a large A77 in Fig. [2] means that one of the 
particles is probably at large r\ in Fig. [5j this evidently is 
not a significant effect on the A77 dependence. 



B. pt dependence 



IV. RESULTS 

First we will show v 3 vs. 77 using two standard event 
plane methods, followed by V3 vs. px for these methods 
and also for the wide Gaussian two-particle correlation. 
Finally, we present the integrated V3 vs. centrality for 
these methods and also for the two-particle cumulant 
method with an 77 gap. Results in all the figures are 
presented with only statistical errors unless stated oth- 
erwise. 



The pt dependence is shown in Fig. [6] For the wide 
Gaussian method, Eq. ^ was used together with the 
parameters from Fig. [3] for each px bin. The results for 
the wide Gaussian method with either kind of weighting 
are almost the same as those for the TPC using sub-event 
planes, meaning that for either of these two methods the 
narrow Gaussian does not significantly affect the wide 
Gaussian. However, in Fig. [7] the results with the event 
plane in the FTPCs are considerably lower, presumably 
because of the larger 77 gap to be discussed in Sec. |IVD| 



A. 77 dependence 

Figure [5] shows the r] dependence of V3 using two event 
plane methods. For particles in the TPC using the 77 
sub-event method, V3 is slightly peaked at mid-rapidity. 



C. Centrality dependence 

Figure [S] shows the centrality dependence of V3 ob- 
tained by integrating over px using the observed yields. 
Shown are two-particle cumulants f3{2} with a minimum 
pseudorapidity separation between particles of one unit. 



7 



0.1 
0.05h 




> 



0.1h 

> 

0.05I- 







0.1- 
0.05- 




0-10% Central 
• v 3 {TPC}:STAR 
v 3 {FTPC}:STAR 



20-30% Central 



40-50% Central t 

• ■• 

Ur. 



10-20% Central 

. • • 



30-40% Central 
. « " • tf • * ^ " 



50-60% Central 
* * 



i .** 



1 2 3 4 1 2 3 4 

P T (GeV/c) 

FIG. 7: (Color online) The third harmonic coefficient as a 
function of pr, for different centralities for Au+Au collisions 
at ^/JM = 200 GeV, for tracks in the TPC with -1.0<tj<1.0. 
The event planes are constructed either in the TPC or in the 
FTPCs. 
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Shown also is «3{2} from Eq. Q and Fig. [2] for the wide 
Gaussian using particle pair weighting. Shown also are 
W3{TPC} and 1>3{FTPC} where u 3 is measured relative 
to the third harmonic event plane reconstructed either in 
the TPC sub-events or the FTPCs. Using unit weights 
in Eq. ([I]) slightly lowered the wide Gaussian results for 
very peripheral collisions. For t>3{2} without a Arj cut 
the curve would be a factor of two higher for peripheral 
collisions and off scale. 



Systematic uncertainties have been estimated by vary- 
ing the DCA track cuts and the number of fit points, 
the event cut of vertex z, and the event plane flatten- 
ing method. These uncertainties have been combined in 
quadrature to obtain the systematic uncertainties shown 
in Fig. [8] The correlat ion of the third and second har- 
monic event planes was investigated by (cos 6(^3 — ^2)) 
and within the statistical uncertainties was found to be 
consistent with zero for this data set. This is reason- 
able for this mixed harmonic result since observing the 
correlation between the third and second harmonic event 
planes requires a three particle correlation analysis to fix 
the direction of the first harmonic event plane [5] . 



FIG. 8: (Color online) The third harmonic coefficient as 
a function of centrality from different methods of analysis 
for Au+Au collisions at ^/sjvjv = 200 GeV, integrated for 
0.15<p T <2.0 GeV/c and -1.0<r?<1.0. The curves connect 
the points and the bands show the systematic uncertainties. 
The systematic errors of the wide Gaussian method are simi- 
lar to those for the TPC event plane method. 



D. A77 dependence 

Clearly the various analysis methods for v 3 differ 
greatly in Fig. [8] The results from the wide Gaussian and 
the TPC event plane are similar, showing that the nar- 
row Gaussian effect is eliminated in both. When a large 
Ar] is specified the v 3 values decrease, especially for the 
peripheral collisions in Fig. [8] The variation between the 
sets of results in Fig. [8] is caused by the A77 dependence 
as shown in Fig. [9] Two-particle correlation results in 
the TPC as a function of A77 for three charge combina- 
tions and two centralities are shown. Also shown are the 
results for three analysis methods as a function of the 
mean Arj of the particles. For the points at |A?j| = 3.21 
the event plane resolutions may be a bit low, and thus 
the V3 values slightly high, because the -q gap between the 
two FTPCs is larger than that between the particles and 
the event plane. There is general agreement in the grad- 
ual decrease of 1^ with Ar). The non-flow contributions 
due to short range correlations are effectively suppressed 
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FIG. 9: (Color online) The square of the third harmonic coefficient as a function of pseudorapidity separation for Au+Au 
collisions at ^snn ~ 200 GeV for tracks with 0.15<pr<2.0 GeV/c. Shown are Unlike Sign, Charge Independent, and Like 
Sign results at 0-5% centrality (open symbols) and 30-40% centrality (closed symbols). Most of the points at low \Ar]\ are not 
plotted because they correspond to the narrow Gaussian and go off the top of the scale. Also shown, by larger symbols, are the 
squares of the mean values (connected by purple lines) at the same two centralities from three analysis methods: The point 
at |Ar/| = 0.63 is from the sub-event method using the TPC with |r/|<l. The point at 1.33 is from the 2-particle cumulant 
method with |Atj|>1. The point at 3.21 is from correlations with the FTPC event plane. The dashed (green) curves without 



points are from a minimum bias Glasma calculation [35 for , 
acceptance with overall normalization set to the data at jAr/| = 



= 200 GeV with 0.15<p T <2.0 GeV/c done for the STAR 



by using either the wide Gaussian or by an rj gap. This 
result is consistent with previous studies of elliptic flow 
based on two-particle correlations, but in a previous work 
the corresponding wide Gaussian was ascribed to mini- 
jet correlations [30]. The decrease with A 77 of f^{2} has 
been seen previously by the ATLAS Collaboration [31]. 
It has been calculated in Ref. [32] as a decrease in non- 
flow. The decreasing effect of fluctuations from initial 
state gluon correlations has been described in Ref. [28] 
but without evolution to the final state. The dilution of 
fluctuations during transport to the final state has been 
calculated in Ref. [33]. Reference [33] also describes the 
decorrelation of flow with increasing pseudorapidity gap 
using the AMPT model. Figure [9] also is reminiscent of 
the well known near-side ridge in a plot of Ar] vs. A<p 
having a peak and shoulder |17j . The far-side ridge may 
also contribute to this shoulder. 

As Fig. [9] shows, we did not find that v 3 stabilized 
at a constant value for large A77 within the acceptance 
of STAR. Thus one might ask if one should extrapolate 
to large Ar] to avoid non-flow, or small Arj to measure 
all the fluctuations. However, it is clear that one must 
always quote A77 for each v% measurement and one must 



compare results to models with approximately the same 
Ar] as the experiment. To help clarify the physics we 
compared like and unlike charge-sign combinations, since 
they have different contributions from resonance decays, 
fluctuations, and final state interactions, but found little 
difference. One source of fluctuations is calculated in 
the Glasma model [35] and shown by the Glasma lines, 
normalized to fit the data at Ar] =1 in the figure. They 
show some decrease with A77, but not as much as in the 
data. 



E. Four-particle cumulants 

The results from four-particle cumulants, «3{4}, with 
weighting by the number of combinations are shown in 
(top). They are consistent with zero within the 



Fig. 10 



errors, in contrast to the ALICE results [35] at the higher 
beam energy. Four-particle cumulants are known to sup- 
press non-flow and Gaussian fluctuations [371 121] • To 
look for non-Gaussian fluctuations, Ref. [39] suggests 
plotting (2 * vl{2} - vl{i})/vl{2} = 2 - «£{4}/«£{2}. 



This ratio, which is shown in Fig. 10 (bottom), on the 



9 



0.15 



» o.ih 
o 



o.o5h 



> 



i * * 

-n 05 - , 



2.4- 



CM 

t CO 
> 



2.2- 



% n 2H 



CM 



1.8- 



• ALICE 
■ STAR 



i 



10 20 30 40 50 
Central ity (%) 



FIG. 10: (Color online) Top: The fourth power of the 
third harmonic coefficient from four-particle cumulants as 
a function of centrality for Au+Au collisions at ^/snn = 
200 GeV, with track selections 0.15<p T <2.0 GeV/c and 
— 1.0<fJ<1.0. The ALICE results [36] are for Pb+Pb colli- 
sions at ^snn = 2.76 TeV, with track selections 0.2<pt<5.0 
GeV /c and — 0.8<r?<0.8. Bottom: The points in the top fig- 
ure divided by the fourth power of the third harmonic flow 
from the r\ sub-event method, showing the deviation from 2. 



event plane the average \Ar]\ was 3.21. The PHENIX 
results used the event plane from their RXN detector 
at an intermediate rj of 1.0<?7<2.8. Our results with 
the event plane in the TPC are very similar to those 
of PHENIX. This is surprising because the mean rj of 
their RXN detector is larger than that for the sub-events 
in our TPC. Our FTPC results in Fig. [7j however, are 
lower than theirs. This is reasonable because the mean 
| Arj\ is considerably larger in the FTPC than in the RXN 
detector. 

Comparison to LHC results for Pb+Pb at ^/i/wv = 
2.76 TeV for ALICE [5S] and ATLAS [3T] are also shown 
in Fig. [IT] ALICE results are for |?7|<0.8 and |A?j|>1.0. 
ATLAS results are for |r/|<2.5 with the event plane in 
the forward calorimeter at 3.2<?y<4.9, giving |A?7|>0.8. 
Agreement is good not only between RHIC experiments, 
but also between RHIC and LHC experiments. This is 
surprising because of the somewhat different At? ranges. 



VI. MODEL COMPARISONS 

In the event-by-event ideal hydro model, V3 was stud- 
ied first by Ref. [IT], and then by Ref. gTJ. Refer- 
ences [321 S3] concluded that instead of averaged initial 
conditions, event-by-event calculations are necessary to 
compare with experimental data. The first prediction of 
7j 3 with viscous hydro was in Ref. [10] . The linear trans- 
lation from initial space fluctuations to final momentum 
fluctuations has been calculated for elliptic flow with the 
NeXSPheRIO model [44]. Reference [22] calculates the 
additional fluctuations induced during the viscous expan- 
sion. 



A. Pseudorapidity separation 



average overlaps with both the ALICE results and the 
Gaussian value of 2. Even though the differential vz{pt) 
values for STAR and ALICE (which will be shown later) 
are the same, the integrated results for ALICE are larger, 
making their error bars in this figure smaller. Also, 
ALICE results come from a higher multiplicity at their 
higher beam energy, probably making the non-Gaussian 
effect more visible. Alternatively, the non-Gaussian fluc- 
tuations only may appear at the higher pj< values in- 
cluded in the ALICE results. However, the precision of 
the STAR data does not allow us to conclude whether 
the STAR fluctuations are Gaussian or not. 



COMPARISONS TO OTHER EXPERIMENTS 



Figure 11 compares our ^{TPC} results from Fig. [7] 
with those from PHENIX [40]. The PHENIX results are 
shown for \r]\ < 0.35, while for STAR the 77 acceptance 
was \r]\ < 1.0. For the STAR results from the TPC the 
mean | A77I was 0.63, while for the results using the FTPC 



Calculations of vf {2} vs. Arj have been done in 
Ref. |32| . They used an event-by-event viscous hydro 
model and addressed the effect of radial flow on local 
charge conservation in hadronization. Their results have 
a similar slope vs. Aij as the data in Fig. [9} but the values 
are too high. The normalization to fit the data probably 
could be adjusted. But their charge balancing mechanism 
would predict a much bigger difference between unlike- 
sign pairs and like-sign pairs. There is only a small spread 
in the data in Fig. [9] at Ai] about 0.5, largely ruling out 
this mechanism. 

The Glasma model calculations of Ref. [33] show some 
decrease in i>| with (A^] in Fig. giving a partial expla- 
nation for the decrease with \Arj\. However, these calcu- 
lations for the initial state are not sufficient to explain the 
sharper fall off of v\ vs. | A77I seen in the data. This perbu- 
tative model is strictly only valid at the higher pt values 
(Pt 3> Qs, where Qs is the saturation scale of the bulk 
matter produced in the collision). Reference [35] says 
"The decorrelation of the two-particle correlation with 
increasing rapidity gap demonstrates the violation of the 
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FIG. 11: (Color online) The third harmonic coefficient as a function of transverse momentum, for different centralities. The 
STAR « 3 {TPC} results are from Fig.[7j Also shown are PHENIX results, ATLAS results starting at 10% centrality, and ALICE 
results for 30-40% centrality. 



boost invariance of the classical Glasma flux tube picture 
by quantum evolution effects." In principle the normal- 
ization could be determined by hydrodynamic transport 
to the final state. However, it is probable that the large 
discrepancy between the methods in Fig. [8] has its origin 
in the Ary dependence of fluctuations, either in the initial 
state or in the hydrodynamic evolution. 

Another Glasma flux tube model with radial flow has 
been used to calculate fluctuations and V3 [45]. Refer- 
ence [T7J says that the near-side ridge caused by long- 
range 77 correlations, and odd harmonics in the azimuthal 
anisotropy, are two ways of describing the same phe- 
nomenon - the response of the system to fluctuations 
in the initial density distribution. 



B. Transverse momentum dependence 



In Fig. 12 



V2 [461 and V3 obtained with the TPC sub- 
event plane method are compared as a function of trans- 
verse momentum with several models for 0-5%, 20-30%, 
and 30-40% central collisions. The experimental results 
for the TPC sub-event plane method are shown because 



they eliminate the short-range correlations but yet have a 
small IA77I like the theory calculations. Shown in Fig. 12 
are the ideal and viscous event-by-event hydrodynamic 
model of Refs. [131 HH] where the initial conditions come 
from a Monte Carlo Glauber model and the ratio of shear 
viscosity (77) to entropy density (s) is r]/s = 0.0, 0.08, 
and 0.16. To properly include fluctuations, 100 to 200 
events were simulated and then the root-mean-square 
flow values calculated. The agreement with the hydro 
for ri/s = 0.08 is very good. NeXSPheRIO [41] root- 
mean-square results for 20-30% and 30-40% centralities 
below one GeV/c are also good. Also shown are the 
results from the AMPT model [H] with string melting 
for the latest set of parameters ("Set B"). The agree- 
ment for i> 2 is good, but the calculated V3 is a bit high. 
AMPT has also been used for v% from symmetric [UJ 
and asymmetric collisions |48) . Predictions for V3 from 
Parton-Hadron-String Dynamics [49] at 30-40% central- 
ity for 1 77 1 <0.5 have been made by the sub-event method 
with the event planes at 1.0<|?7|<4.0, and show good 
agreement in the figure lower right. HIJING [50] does 
not predict any significant V3 as u| in the pt range up to 
1.5 GeV/c is both negative and positive, with absolute 
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FIG. 12: (Color online) v 2 (top) and v 3 (bottom) for Au+Au collisions at = 200 GeV in 0-5% (left), 20-30% (middle), 

and 30-40% (right) centrality as a function of transverse momentum, compared with ideal and viscous hydro [T3], AMPT 
transport [14], NeXSPheRIO |41| . and Parton-Hardon-String Dynamics |49| models. The STAR v 2 values (top) are from 
Ref. 06]. 



values less than 2 x 10~ 4 . 

Previously, elliptic flow results were mostly described 
by hydro with r\j s = 0.08 with Glauber initial conditions 
in the case of mid-central collisions [T3j . We find that i>3 
results are described by this model with a similar viscos- 
ity. The NeXSPheRIO model at low p T and the PHSD 
model also agree with the data. 

VII. SUMMARY 

We have presented STAR's measurements of third har- 
monic flow of charged particles from Au+Au collisions at 
y / ijvjv = 200 GeV as a function of pseudorapidity, trans- 
verse momentum, rapidity gap, charge sign, and central- 
ity. We report results from a two-particle method for 
particle pairs with an ?y gap or fit with a wide Gaus- 
sian in pseudorapidity separation, as well as from the 
standard event-plane method with the event plane near 
mid-rapidity or at forward rapidity. Short-range corre- 
lations are eliminated either by an rj gap or by discard- 
ing the narrow Gaussian in pseudorapidity separation. 
The measured values of V3 continuously decrease as the 
mean pseudorapidity separation of the particles increases 
within the range observable by STAR. A model for non- 
flow predicts a big difference between different charge 
sign pairs which is not observed in the data. A model for 
the decrease of fluctuations with pseudorapidity separa- 
tion from a Glasma initial state has some aspects of the 



data. Because of this, and the good agreement of vz{pr) 
with models including fluctuations, it is likely that V3 is 
mainly due to Ar] dependent fluctuations [2S]. Accord- 
ing to the models these fluctuations should be largely 
independent of beam energy. 
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